arXiv: astro-ph/0512176v2 17 Apr 2006 


Evidence of a Change in the Long Term 
Spin-down Rate of the X-ray Pulsar 4U 1907+09 

A. Baykal 1 , S.Q. inam 2 , E. Beklen 1,3 
1 Physics Department, 

Middle East Technical University, 06531 Ankara, Turkey 
altan@astroa.physics.metu.edu.tr 
2 Department of Electrical and Electronics Engineering, 
Ba§kent University, 06530 Ankara, Turkey 
inam@baskent.edu.tr 
3 Physics Department, 

Suleyman Demirel University, 32260 Isparta, Turkey 
elif@astroa.physics.metu.edu.tr 


Abstract 

We analyzed RXTE archival observations of 4U 1907+09 between 17 
February 1996 and 6 March 2002. The pulse timing analysis showed that 
the source stayed at almost constant period around August 1998 and then 
started to spin-down at a rate of (—1.887 + 0.042) x 10~ 14 Hz s _ l which is 
~ 0.60 times lower than the long term (~ 15 years) spin-down rate (Baykal 
et al. 2001). Our pulse frequency measurements for the first time resolved 
significant spin-down rate variations since the discovery of the source. We 
also presented orbital phase resolved X-ray spectra during two stable spin 
down episodes during November 1996 - December 1997 and March 2001 

- March 2002. The source has been known to have two orbitally locked 
flares. We found that X-ray flux and spectral parameters except Hydrogen 
column density agreed with each other during the flares. We interpreted the 
similar values of X-ray fluxes as an indication of the fact that the source 
accretes not only via transient retrograde accretion disc (in’t Zand et al. 
1998) but also via the stellar wind of the companion (Roberts et al. 2001), 
so that the variation of the accretion rate from the disc does not cause 
significant variation in the observed X-ray flux. Lack of significant change 
in spectral parameters except Hydrogen column density was interpreted as 
a sign of the fact that the change in the spin-down rate of the source was 
not accompanied by a significant variation in the accretion geometry. 
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1 Introduction 

The X-ray source 4U 1907+09 is an accretion powered X-ray pulsar accreting 
plasma from a blue supergiant companion star. Since its discovery in the Uhuru 
survey (Giacconi et al. 1971), it has been observed by various X-ray observatories 
including Ariel V (Marshall & Ricketts 1980), Tenma (Makishima et al. 1984), 
EXOSAT (Cook & Page 1987), Ginga (Makishima & Mihara 1992; Mihara 1995), 
BeppoSAX (Cusumano et al. 1998), XMPC (Chitnis et al. 1993), IXAE (Muk- 
erjee at al.2001) and RXTE (in’t Zand, Baykal & Strohmayer 1998a; in’t Zand, 
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Strohmayer & Baykal 1997, 1998b; Roberts et al. 2001 ;Baykal et al. 2001). 
Marshall & Ricketts (1980) determined the orbital period to be 8.38d using Ariel 
V observations. They also found two flares, namely a primary and a secondary 
occuring near periastron and apastron, and EXOSAT and RXTE observa¬ 
tions (in’t Zand et al. 1998a,b) have clearly shown that these flares are locked 
to orbital phases separated by half an orbital period. Using Tenrna observations, 
Makishima et al. (1984) also found the pulse period of the source to be 437.5 
sec. Makishima et al. (1984) and Cook & Page (1987) discussed the possibility 
that the two flares are caused by an equatorial disc-like envelope (circumstcllar 
disc) around the companion star. When the neutron star crosses the disc, the 
mass accretion rate onto the neutron star increases temporarily. We expect to 
see two peaks while the neutron star passes twice through the circumstellar disc 
of the companion star in every cycle. We also expect a transient accretion disc 
formation around the neutron star while the neutron star passes through this 
circumstellar disc. For the formation of the circumstcllar disc around the com¬ 
panion star, the companion must most likely be a Be type star. Recent optical 

observations showed that the optical companion of 4U 1907+09 is an 08/09 su¬ 
pergiant with a dense stellar wind (Cox, Kaper, Mokiern 2005). Formation of a 
transient accretion disc around the neutron star as a result of disruption of the 
dense stellar wind of a massive companion is quite possible as well (Blondin et 
al. 1990). 

The timing measurements from EXOSAT (Cook & Page 1987) and RXTE 
observations (in’t Zand et al. 1998a,b; Baykal et al. 2001) showed continuous 
spin down of the neutron star.The source was found to be spinning down almost 
at a rate of v — (—3.54 ± 0.02) x 10~ 14 i7+s _1 for more than 15 years (Baykal et 
al. 2001). Timing analysis of RXTE (in’t Zand et al. 1998) and IXAE (Mukerjee 
et al. 2001) observations of 4U 1907-09 also revealed transient oscillations with 
periods of about 18.2s and 14.4s respectively. These transient oscillations were 
suggested to be related to the presence of a transient retrograde accretion disc 
which is also responsible for the slowing down of the pulsar. 

In this paper we present the results of timing and spectral studies of archival 
RXTE-PCA observations of 4U 1907+09, and report a change in the spin down 
rate of this system. 

2 Instrument and Observations 

The Proportional Counter Array (PCA)(Jahoda et al.1996) on board the RXTE 
consists of five identical proportional counter units (PCUs) coaligned to the same 
point in the sky, with a total effective area of approximately 6250 cm 2 , and a held 
of view of ~ 1° FWHM, operating in the 2-60 keV energy range. The number of 
active PCUs is varied between one and five during the observations.Observations 
after 2000 May 13 belong to the observational epoch for which background level 
for one of the PCUs (PCU0) increased due to the fact that this PCU started 
to operate without a propane layer. Latest combined background models (CMs) 
were used together with FTOOLS 6.0 to estimate the appropriate background. 

The RXTE-PCA observations used in this work are listed in Table 1. 
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3 Analysis and Results 

3.1 Timing Analysis 

In order to obtain background subtracted light curves, we generated the back¬ 
ground using the background estimator models based on the rate of very large 
solar events, spacecraft activation and cosmic X-ray emission with the standard 
PCA analysis tools. Then we subtracted background light curves from the source 
light curves obtained from the event data. The background subtracted light 
curves were also corrected to the barycenter of the solar system and to the binary 
orbital motion of 4U 1907+09 around its companion. For the correction of binary 
orbital motion, we used the orbital parameters deduced by In’t Zand, Baykal & 
Strohmayer (1998a). In the timing analysis we fold the light curve outside the 
intensity dips from long data string around 7-10 ksec on statistically independent 
trial periods (Leahy et al. 1983). A template pulse profile was constructed by 
folding the data on the period giving the maximum y 2 . Template pulse profile 
was represented by its Fourier harmonics (Deeter & Boynton 1985) and cross- 
correlated with the harmonic representation of average pulse profiles from each 
observation. The pulse phases obtained from the cross correlation analysis are 
presented in Figure 1. 

In order to estimate pulse frequency derivatives, the pulse phases are fitted to 
the quadratic polynomial 

S(j)= (j) Q + 5v(t - t Q ) + iz>(£ - t 0 ) 2 (1) 

where 5(j) is the pulse phase offset deduced from the pulse timing analysis, t 0 is the 
mid-time of the observation, <f> 0 is the phase offset at t 0 , 5v is the deviation from 
the mean pulse frequency (or additive correction to the pulse frequency), and v is 
the pulse frequency derivative of the source. The pulse phases and the residuals of 
the fit after the removal of the quadratic polynomial are also presented in Figure 
1. Table 2 presents the timing solution of 4U 1907+09. 

The pulse frequency derivative from the sequence of 19 pulse phases spread 
over 360 days was measured as v = ( — 1.887 ± 0.042) x 10 -14 Hz s _1 . This 
value approximately ~ 0.60 times lower then previous measurements of pulse 
frequency derivatives which was measured over 383 days times span (Baykal et 
al. 2001). In this and the previous observations Baykal et ah, (2001) time span 
of the observations are similar (~ over year) and in both observations possess 
low timing noise. However the spin down rate is lowered by a factor of 0.6 times 
in latter observations. 

In order to see the spin down trends clearly we estimate the pulse frequency 
histories by taking the derivatives of each pairs of pulse phases from the pulse 


Table 1: Observation list for 4U 1907+09 


Time of Observation 
start end 

dd mm yyyy dd mm yyyy 

XTE Prop. ID 

IN umber of 
Observations 

17 02 1996-23 02 1996 

10155 

8 

25 11 1996-14 12 1997 

10154-20146 

24 

26 07 1998-01 10 1998 

30093 

14 

10 03 2001-06 03 2002 

60061 

38 
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Figure 1: Pulse phase and its residuals fitted to the orbital model presented in 
Table 2. 


phases of this work and previous work by Baykal et ah, (2001) and presented in 
Table 3 and Figure 2 along with previous frequency measurements. 

The whole dataset of 4U 1907+09 mentioned in previous sections was also 
analyzed for possible detection of any transient oscillations around ~ 14 — 18s, 
similar to those reported from previous RXTE (in’t Zand et al. 1998a) and IXAE 
(Mukerjee at al.2001) observations. In the power spectra of 3-25 keV RXTE- 
PCA light curve, we searched for transient oscillations around ~ 0.06 Hz in the 
0.50s binned ~ 170 light curve segments of length of 512s. In order to test the 
significance of these oscillations, we used the approach of van der Klis (1989). As 
the result of this search, we found only 3 512s segments in the whole dataset with 
the transient oscillations of significance exceeding 5cr at MJD ~ 50136, ~ 51079 
and ~ 51087. Our results about transient oscillations are marginal, and therefore 
inconclusive. Future observations may be useful to observe and analyze possible 
transient oscillations of the source. 

3.2 Analysis of the Energy Spectrum 

The complete dataset of 4U 1907+09 mentioned in Section 2.1 were also analyzed 
to investigate whether there is any significant spectral evolution accompanied 
with the change in the spin down rate of the pulsar or not. Spectra, background 
and response matrix hies were created using FTOOLS 6.0 data analysis software. 
Background spectra were generated using the background estimator models based 
on the rate of very large events, spacecraft activation and cosmic X-ray emission. 
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Table 2: Timing solution of 4U 1907+09 a 


T n /2 (Orbital Epoch) (MJD) 

50134.76 + 0.06 b 

Porb(d) 

8.3753 + 0.0001 b 

o^sin i (It s) 

83 + 2 b 

e 

0.28 + 0.04 b 

w 

330 + 7 6 

t 0 (MJD) c 

52154.6183 + 0.0005 

Pulse Period (s) 

441.0932 + 0.0003 

Pulse Period derivative (s s“l) 

(3.653 + 0.081) x 10“ 9 

Pulse freq. derivative (Hz s _ l) 

(-1.887 + 0.042) x 10“ 14 


a Confidence intervals are quoted at la level. 

b Orbital parameters are taken from In’t Zand et al. (1998a).P Q r&—orbital 
period,a^sin i=projected semimajor axis, e=eccentricity, w=longitude of peri- 
astron. 

c Mid-time of the observation used in pulse timing (see Equation 1 in 
the text) 


Spectral analysis and error estimation of the spectral parameters were per¬ 
formed using XSPEC version 11.3.2. We extracted individual 3-25 keV PCA 
spectra each corresponding to a single pulse period. For these spectra, energies 
lower than 3 keV were ignored due to uncertainties in background modelling while 
energies higher than 25 keV were ignored as a result of poor counting statistics. In 
order to take into account the X-ray background from the Galactic ridge and the 
nearby supernova remnant W49B, we followed the approach used by Roberts et al. 
(2001): We extracted overall ’’dip” state spectra of the source for each of the pro¬ 
posal IDs listed in Table 1 and the corresponding background spectra. We found 
that the count rate and the spectral shape of the dip state spectra were consistent 
with the models of the diffuse emission from the Galactic ridge (Valinia & Mar¬ 
shall 1998). So, these dip state observations were used as additional background 
spectra in spectral fits of the source spectrum. These additional background spec¬ 
tra were created by subtracting the corresponding dip state model background 
spectrum from the dip state source spectrum using ’’mathpha” FTOOL. 

A power-law model with low-energy absorption (Morrison & McCammon 
1983), multiplied by an exponential high energy cut-off function (White, Swank 
& Holt 1983) was used to model 3-25 keV spectra of the source. We found no 
evidence of a statistically significant iron line emission at ~ 6.4 keV which is con¬ 
sistent with the previous RXTE spectral fits (Roberts et al. 2001). An additional 
cyclotron resonance absorption line (Makishima et al. 1990) component at ~ 19 
keV was required to fit the spectra. The parameters of cyclotron absorption line 
was found to be consistent with the fundamental cyclotron line parameters found 
from BeppoSAX observations of the source (Casumano et al. 1998). Obtaining 
the spectral fits, we fixed the cyclotron absorption line energy at 19 keV. The 
result of the spectral fits as a function of orbital phase were presented in Figure 
3. Sample model parameters were listed in Table 4. 
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Figure 2: Plot of pulse period history listed in Table 3 (Makishima et al. 1984, 
Cook & Page 1987, Mihara 1995, In’t Zand et al. 1998, Baykal et al. 2001 and 
this work). Solid line and the residuals correspond to the previous spin-down rate 
found by Baykal et al. (2001). 


4 Conclusion 

The spin-down rate of 4U 1907+09 obtained from the pulse arrival times between 
March 2001 and March 2002 was found to be ~ 0.60 times lower than both 
the previous RXTE measurements of spin-down rate between November 1996 - 
December 1997 by Baykal et al., (2001) and the long term spin down trend of the 
source between 1983 and 1997 (See Figure 1 in Baykal et al. 2001). In this and 
the previous RXTE observations, time span of the observations were similar (~ 
over year) and in both observations posses low timing noise. However spin down 
rate is lowered by a factor of 0.6 in latter observations. It is also interesting to 
observe that the spin-down rate is consistent with zero around MJD 51000 before 
this significant change in spin-down rate. Our pulse frequency measurements for 
the first time resolved significant spin-down rate variations since the discovery of 
the source. 

The steady spin-down rate (Baykal et al. 2001) and the presence of occasional 
transient oscillations (in’t Zand et al. 1998a; Mukerjee et al. 2001) in 4U 1907+09 
support the idea that the source accretes from retrograde transient accretion disc. 
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The possibility of a prograde disc for which the magnetospheric radius should be 
close to the corotation radius so that the magnetic torque overcomes material 
torqne is not likely to be the case for this system (in’t Zand et al. 1998a). If 
the disc is retrograde and the material torqne is dominant, using Ghosh& Lamb 
(1979) accretion disc model, a decrease in spin-down rate of the neutron star 
should be a sign of a decrease in the mass accretion rate coming from the accretion 
disc. If the disc accretion is the only accretion mechanism, this decrease in mass 
accretion rate should also lead to a decrease in X-ray flux of the system. 

From the X-ray spectral analysis of the source, we found no clear evidence of 
a correlation between 3-25 keV flux and the change in spin-down rate (see Figure 
3). The flux levels were about the same for latter observations with low spin- 
down rate compared to the rest of the observations. This was not our expectation 
if we only consider disc accretion and may be because of the fact that the total 
mass accretion rate is not only due to only disc accretion, but also accretion from 
stellar wind may contribute to the total mass accretion rate (see also Roberts et 
al. 2001) . In case of accretion from both transient disc and wind, the change in 
mass accretion rate from the accretion disc might not cause a substantial change 
in X-ray flux. However, it is important to note that the orbital coverage of the 
latter (proposal ID 60061) observations is poor and limited to the phase locked 
flares in orbital phases ~ 0.05 and ~ 0.6 as seen in Figure 3. Future observations 
of the source may be helpful to have a better understanding of a possible relation 
between spin-down rate and X-ray flux. 

Examining spectral parameters plotted in Figure 3 especially for flare parts, 
we also found an evidence of an increase in Hydrogen column density (nn) for the 
latter observation with lower spin-down rate. On the other hand, there is no clear 
evidence of a variation of the other spectral parameters in the flares after the spin- 
down rate decreased. Change in nn may be related to a change in the accretion 
geometry of the source. However, if there had been an accretion geometry change 
accompanied with the change in the spin-down rate of the source, we would expect 
other spectral parameters to vary as well, especially the power law index as in 
the RXTE observations of SAX J2103.5+4545 (Baykal et al. 2002), and 2S 1417- 
62 (Inam et al. 2004). New X-ray observations of the source will be useful to 
monitor any possible changes in the spin-down rate and its consequences. 
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Table 3: RXTE pulse period measurements of 4U 1907+09 


Epoch Pulse Period 
(MJD) (s) 

45576 437.483+0.004 

45850 437.649+0.019 

48156.6 439.19+0.02 

50134 440.341+0.014 

50424.3 440.4854+0.0109 

50440.4 440.4877+0.0085 

50460.9 440.5116+0.0075 

50502.1 440.5518+0.0053 

50547.1 440.5681+0.0064 

50581.1 440.5794+0.0097 
50606.0 440.6003+0.0115 

50631.9 440.6189+0.0089 

50665.5 440.6323+0.0069 

50699.4 440.6460+0.0087 

50726.8 440.6595+0.0105 

50754.1 440.6785+0.0088 

50782.5 440.6910+0.0097 

51021.9 440.7045+0.0032 

51080.9 440.7598+0.0010 

51993.8 441.0484+0.0072 

52016.8 441.0583+0.0071 

52061.5 441.0595+0.0063 
52088.0 441.0650+0.0063 

52117.4 441.0821+0.0062 

52141.2 441.0853+0.0082 

52191.4 441.1067+0.0046 

52217.2 441.1072+0.0077 

52254.3 441.1259+0.0074 
52292.0 441.1468+0.0065 

52328.8 441.1353+0.0090 


Reference 


Makishima et al.1984 
Cook & Page 1987 
Mihara 1995 
In’t Zand et al.1998 
This work 
Baykal et al.2001 
This work 
This work 
This work 
This work 
This work 
This work 
This work 
This work 
This work 
This work 
This work 
Baykal et al.2001 
Baykal et al.2001 
This work 
This work 
This work 
This work 
This work 
This work 
This work 
This work 
This work 
This work 
This work 


Table 4: Sample Spectral Parameters of Individual PCA Observations of 4U 
1907+09 


Parameter 

iVIJD 50547 

MJD 52061 

Hydrogen Column Density (10 22 cm -2 ) 
Power Law Photon Index 

Power Law Norm. (10 _2 .cts.cm _2 .s _1 ) 
Cut-off Energy (keV) 

E-folding Energy (keV) 
Cyclotron Depth 

Cyclotron Energy (keV) 
Cyclotron Width (keV) 

3-25 keV LInabsorbed X-ray Flux 
(10~ lo ergs.cm~ 2 .s -1 ) 

Reduced y 2 (43 d.o.f.) 

5.51 =F 0.50 
1.20 =F 0.04 
4.04 =|= 0.09 
12.9 =F 0.2 
7.0 =F 0.5 
0.14 =|= 0.05 
19.0 (fixed) 

1.12 =f 0.10 

6.36 =F 0.14 

0.74 

9.13 =F 0.23 
1.21 =|= 0.03 
2.83 =F 0.15 
13.7 =F 0.3 
9.6 =F 0.5 
0.20 =F 0.02 
19.0 (fixed) 
1.71 =F 0.34 
4.69 =F 0.25 

0.98 
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Figure 3: Evolution of 3-25 keV unabsorbed flux, Hydrogen column density, 
power law index, cut-off energy, e-folding energy, and reduced y 2 in orbital phase. 
Vertical dashed lines indicate the orbital phase corresponding to periastron pas¬ 
sages. Points marked with triangles indicate observations between March 2001 
and March 2002 (proposal ID 60061) for which spin-down rate of the source was 
lower. The dotted lines indicate ±er from the periastron. 
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